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Abstract 

X^/^ \ We calculate the distribution of the fc-th largest eigenvalue in the random matrix Levi-Smirnov 

(LSE) ensemble, using the spectral dualism between LSE and chiral Gaussian Unitary Ensemble 
(GUE). Then we reconstruct universal spectral oscillations and we investigate an asymptotic behavior 
of the spectral distribution. 

^ ; 1 Introduction 
O 

Q\ ' Random matrix ensembles provide a powerful and generic formalism allowing addressing several statistical 

problems of the energy spectra in various complex quantum systems. In most applications, the ensembles 
used belong to the basin of attraction of the Gaussian ensembles. Relatively few is known about the 
mathematical properties of the matrix ensembles with so strongly fluctuating elements, that the average 
and the variance of the ensemble diverge. Such ensembles can be viewed as the matrix generalization 
of the stable probability distributions of the probability theory, usually known as Levy distributions. 
Historically, such ensembles of matrices were proposed by Bouchaud and Cizeau jjj. Matrix elements 
of the corresponding matrices were drawn from stable, one-dimensional probability distributions. This 
construction, however, was breaking the rotational invariance of the ensemble, imposed in most Gaussian- 
like and polynomial matrix models. Recently, borrowing from the mathematical concepts of free random 
variables analysis ||], a new class of stable, rotationally invariant matrix models was introduced in ||. 
Exploiting the Coulomb gas analogy, the authors ||] propose a general method of constructing explicit 
matrix ensembles characterized by stable power spectra, with asymptotic behavior A~ 1_Q . The stability 
condition means, that the matrix convolution of the identical independent ensembles with power-like 
spectra exhibits, modulo the rescaling and shift, the same power-like spectral distribution. In this sense, 
the ensembles provide a generalization of central limit theorem of Gaussian ensembles, where e.g. the 
convolution of two independent identical Gaussian Unitary Ensembles (GUE) with unit variance is also 
the Gaussian Unitary Ensemble, however with the semi-circle law rescaled by y/2 in the longitudinal 
direction. Due to the growing interest in stochastic processes with long (sometimes called also heavy 
or fat) tails, intermittency and anomalous diffusion, the theory of stable random matrices may form an 
important role in describing generic features of complex systems using new concepts of universality. Very 
recently, new multi-critical random ensembles labeled by continuous scaling exponents were proposed || , 
with eigenvalue density near zero behaving like \x\ a . Similar continuous series of universality class for 
eigenvalues appear naturally in stable random matrix ensembles with a non-compact support It is 
plausible to conjecture strongly [Q, that by simple mapping the scaling behavior of small eigenvalues of 
the continuous multi-fractal regime can be related to the large eigenvalue behavior of the Levy random 
matrix models. 
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In this work, we study one of the simplest ensembles of Levy matrices, so-called Levy-Smirnov matrix 
model. The name origins from the fact, that the model constitutes a matrix analog of the Levy-Smirnov 
probability distribution (j^j), i.e. the spectrum for large eigenvalues is power-like with a — 1/2. In 
the next chapter, we remind the definition and some properties of the LS ensemble. Then, we exploit 
the duality of this model ||] to the massless chiral Gaussian Unitary Ensemble Since the chiral 
GUE plays a crucial role in describing the universal properties of low-lying eigenvalues of QCD Dirac 
operator, an impressive number of analytical tools and methods was developed in few last years. Using 
the aforementioned duality and some of these methods [pi , we calculate the distribution of k-th largest 
eigenvalue in the Levy-Smirnov ensemble (hereafter LSE), and show explicit results for k = 1,2,3,4. 
Finally, we analyze an asymptotic behavior of the spectral function, demonstrating how this behavior 
can be obtained from the microscopic spectral function. We close the paper with the discussion. 



2 Levi-Smirnov ensemble. 

The Random Matrix Theory is usually defined by the partition function 

' d Me- NTlV ^ M \ (1) 



where M is N by N matrix and V(M) is a potential, which, in general, can be even a non-analytical 
function of M. Spectral distribution of the ensemble is defined as 

p(A) = 1 < TtS(X - M) > (2) 

where the averaging is done with the measure defined by the partition function. The Gaussian Unitary 
Ensemble (GUE) is characterized by quadratic potential V(M) = 1/2M 2 (we put variance to 1), leading 
to Wigner semi-circle law for the eigenvalues, p(A) = l/27r\/4 — A 2 , i.e. the eigenvalues are localized on 
the compact interval [—2,2]. 

Levy-Smirnov ensemble is characterized by a rather non-trivial potential || 

V(M) = exp[-i\TTr(M -1 + In M)] (3) 

leading to the spectral function 

Pls(X) = (4) 

localized on the non-compact support [1/4, oo], with asymptotic behavior for large spectra 1/A 1+Q , with 
a = 1/2. This behavior explains the name of the ensemble, since the probability distribution function 
for stable density with similar a and similar maximal asymmetry bears this name 

p L s(x) = -^exp(-l/(2a;))x- 3 / 2 x > (5) 

\/27r 



Note that LS matrix ensemble is not obtained simply by populating the matrix elements from the prob- 
ability distribution (jjj). 



3 A duality between the Chiral and LS ensembles. 

The LS ensemble is one of the simplest example of Levi matrices, because there exist an exact, analytic 
expression for V(M). This model has an unitary symmetry, so we can diagonalize M — > UAU' and, after 
integration over the angular variables, rephrase the partition in terms of the eigenvalues only, the last 
term being the Vandermonde determinant coming from the Jacobian ( see e.g. Q for standard discussion) 

f°° / e -N/\i\ 

z ^ / n^i-rpr- n^-^) • ( 6 ) 
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In this section we show that the distribution of eigenvalues in LSE behaves, up to trivial transfor- 
mation, in the same way as a distribution of the eigenvalues for the chiral Gaussian Unitary ensemble. 
We will show, that the equations describing both LS and Gaussian ensembles can be proved by the same 
method, using almost identical expressions. 

We assume, that N eigenvalues are ordered (Ai > A2 > . . . > Xn)- An unnormalized joint probability 
density f2jv,fc defined as in H in our case takes a form: 

1 r*k N / p -N/\i\ 

i — 1 * i ^> J 

where J\f stands for the normalization: 

1 r°° N / P -«Ai\ 

^ i— 1 \ i / iyj 

The main difference between eq. (Q) and its analog in chGUE lies in the range of integration. In LS case 
we are interested by the large eigenvalues so we integrate over an area with eigenvalues smaller than Afc. 

Using methods originally developed for chGUE in we can calculate the unnormalized joint prob- 
ability taking a microscopic limit. To do this, we change variables in the integral: Ai = l/xi, for i = 
k + 1, . . . , N and we shift variables Xi — ► Xi + 1/A&. This leads us to the expression: 



^« a *» = 7]v3Tv e "^n^2]^n( A i- A i) • (9) 

v >' 1=1 i i<j 



N fc-1 



■/ I] dxie-^XlxlUil-XjXi-Xj/Xk) 2 J] {x t ~x f. 

J ° i=k+l 3 = 1 i>j>k+l 

Now we take a microscopic limit. Motivated by |?], t| [| we go with both Ai and N to infinity keeping 
& = Xi/AN 2 fixed. It is useful to change notation to Ai = l/£i, which help us to find the analogy 
mentioned before. We calculate rescaled function 0,% called (according to M) Wfc(Ai, . . . , Afc): 

/ 4JV 2 \ „ /47V 2 4iV 2 47V 2 



wfc(Ai, . . . , Afc) = ^^J—J^^—,— ...,_j= (10) 

(n ^) -^T^i n^) 2 n ( a * - a2 ) 2 ■ ^-^ + ^ x ( id 
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,00 W / 2 \ 2 w 

x/ n ^^nh-zfe n (^-^ 

" / ° i=fc+i j=i V / i>j>/t+i 

where we have introduced /i 2 = A 2 — A 2 . 

Now, we are able to write this expression using the partition functions for chGUE Z in the microscopic 
limit: 

Uh (A u . . . , A*) = const e^ 4 n A, TT (A? - A 2 ) 2 ^■/'^^■.•■.W-i.W-i) (12) 

2—1 1<J v ' 

This expression shows the deep analogy between chGUE and LS ensembles. With the help of the above 
expression, in the next section, we will derive the rescaled joint probability distribution of the k largest 
eigenvalues in the LSE. 
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4 Distribution of the k-th eigenvalue 

In this section we derive an analytical expression for the distribution of fc-th eigenvalue. This distribution 
is defined as 

pA pA pA 

p k (A)= dAi dA 2 ... dAfc_i w fc (Ai,...,A fc _i,A). (13) 

JO JAi 



Using ( |l2|) and changing variables Xi = \J A\ — A 2 k for i = 1, ... fc — 1 we can write the distribution in 
the form: 



/■A j-Xl rX k - 2 2 k k ~ 1 

Pfe(A) = const / dx\ I dx2 ■ • • / dxk-i e _A / 4 TT Xi (a 

Jo Jo Jo „•_-, 



(z 2 - x 2 ) Z 2 (xi,xi, . . .,x k -i,Xk-x). 

2—1 i<j 

(14) 

An expression showed above contains the partition function for massive chiral GUE model, with 
twofold degeneracy for masses. With the help of the formula for this partition function, proven in we 
can rewrite our final result. In our variables, equation (2.1) from |^| reads 

fc-i 



TT ((a; 2 - a; 2 ) Zf k '(xi, x\, . . . , x k -i, Xk-i) — const det z!f\x a ,x b ) . (15) 

- 1 - J l<a,b<k L 

Using ( |l5| ) we write down our final result: 

/■A rti r%k-2 2 k V 1 

Pfe(A) = const / dx\ I dx 2 ■ ■ ■ / dx k -i e~ A ^ 4 TT xi det Z2 (x ,a;(,) . (16) 

J0 Jo Jo fj{ l<a b< k L J 

where the partition function for two masses was already calculated in 

~(2), v _ x b I 2 (x a )I 3 (x b ) - x a I 2 (x b )I 3 (x a ) 

^2 [x a ,x b )- + ^ ■ \ L <) 

Now, we can go back to the original rescaled variable £ (proportional to un-scaled A). Using substi- 
tution A = -2= we can write the expression for the probability density of the fc-th eigenvalue: 

const [ 2 /^ f xi r^k-i k ~l r 1 

Pk(0 = —&- dx i dx 2 ... dx k -i e^ 1 ^ TT Xi det Z\ \x a , x b ) . (18) 

£ Jo Jo Jo f=l ^ a ' b ^ k 1 J 

This probability needs to be normalized. We calculate this expression for several cases fc = 1,2,3,4. 
The case fc = 1 is very simple: 

Pi(0 = ±e-Vt. (19) 

For higher fc's (greater than 2) calculations are more complicated and cannot be performed by analytical 
methods, so we use the numerical method instead. 

Figure 1 shows the result for probability distribution for first four largest eigenvalues of the LS 
ensemble. 

By definition, the sum of p k distributions reproduces the microscopic spectral density for LSE: 



ps(o = f>(0 = (o- 2 (joivVi) + J?(vVi)) ■ (20) 
fc=i 

calculated in Q. Note also, that the change of variables rewrites this density in terms of microscopic 
spectral density of chiral GUE, calculated by 1 1 



As an independent check of our formulae, we can approximately reconstruct p(£) by adding first four 
distributions. Because the p(£) asymptotic behavior is like l/£ 2 , it is more convenient to plot p-t; 2 instead 
p. Figure 2 shows how oscillatory universal pattern for large eigenvalues is reproduced by the sum of the 
four largest eigenvalues. 
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5 Asymptotic behavior 

In this section, for completeness, we reconstruct the asymptotic behavior of the distribution function p(X) 
for LS ensemble, given by (Q). The asymptotic behavior in infinity is A~ 3 / 2 : 

PLs{X) = 2^^T~ w ^ • ( 21 ) 

We can reconstruct this result from equation (^0|), which we got after taking a microscopic limit. This 
means that we blow up the region close to infinity, because we took the limit A — ► oo keeping fixed. 
This corresponds to asymptotic behavior of the scaling functions for £ close to 0. We use well-known 
property of the Bessel functions (for x real) 

JJx) fa J— (cos(x - -VTT- -) + 0(xr 1 ) ) , for x « 0. (22) 
V 7ra; \ 2 4 / 

Applying this result to the ( ppf ) we obtain the expected result: 

1 Jl ( 2 7T, 2/ 2 7T \ 1 1 



6 Discussion 



The main result of this paper is the equation ( [lq ) , which allows the calculations of pk (£) in an explicit way. 
The analytical form of these probability distributions can be used for comparison to similar empirical 
probability distributions of the largest eigenvalues of stable matrix ensembles, alike it is done in QCD 
lattice analysis for the smallest eigenvalues [ prf . To obtain this result we used well-known methods in 
chiral random matrix models 0, q| and we exploited a non-trivial duality between the chiral Gaussian 
ensemble and LS ensemble. Naively, both ensembles are very different in nature. Chiral Gaussian 
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ensemble is governed by a simple quadratic potential and exhibits a non-trivial symmetry due to the 
off-diagonal block structure. Its spectrum is localized on the compact support, with novel (comparing to 
non-chiral Gaussian ensemble) microscopic oscillations in the vicinity of zero. LS ensemble belongs to non- 
Gaussian stability region, potential is a non- analytic function of M, and spectrum has an infinite support. 
On the basis of the duality known between a and 1 /a distributions for stable one-dimensional probability 
distributions [fl2f one may expect the similar duality between the bulk spectral function for GUE and 
LSE. Indeed, the semi-circle law transmutes into (Q), by changing A 2 into 1/A. This transformation misses 
however the chiral character of the dual ensemble. It is the equivalence of the partition functions at the 
microscopic level, that allows to transmute all non-trivial spectral properties of the chiral ensemble in the 
vicinity of zero to the spectral, universal properties of LS ensemble for large eigenvalues. In some sense, 
LS model is unique, allowing, despite having divergent mean and variance, to use the methods reserved 
for finite variance and compact support models. In general, the study of microscopic properties of Levy 
ensembles requires mapping described in and probably gives only approximate universal scalings. It 
is interesting to study other spectral statistics of the Levy ensembles, especially in connection with the 
critical level statistics in disordered media [jl3| . 
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